GeV-TeV gamma-rays and neutrinos from the Nova V407 Cygni 
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The Fermi-LAT telescope has unexpectedly discovered GeV 7-ray emission from the symbiotic 
Nova V407 Cygni. We investigate the radiation processes due to electrons and hadrons accelerated 
during the explosion of this Nova. We consider a scenario in which GeV 7-ray emission observed 
by Fermi is produced by the electrons with energies of a few tens of GeV in the inverse Compton 
scattering of stellar radiation. On the other hand, the hadrons are expected to reach larger energies, 
due to the lack of radiation losses during acceleration process, producing TeV 7-rays and neutrinos. 
We predict the fluxes of very high energy 7-rays and neutrinos from Novae of the V407 type for 
two models of hadron acceleration and discuss their possible detectability by the present and future 
telescopes (e.g. IceCube, CTA). 

PACS numbers: 97.80. Jp,98.70.Rz, 97.30.Qt, 97.80. Gm 



I. INTRODUCTION 

Nova V407 Cygni belongs to the class of symbiotic No- 
vae which contains a white dwarf and a red giant (RG) 
star of the Mira-type. The binary period of this system is 
43 yrs [l[, and the stars are separated by ~10 14 cm. RG 
has a radius ^500 Hq, the mass loss rate 3 x 10~ 7 M 
yr , the wind velocity ~10 km s _1 and the luminosity of 
10 4 L Q The matter of the RG wind is accreted onto 
the surface of a white dwarf. It is commonly accepted 
that Nova outbursts occur on the surface of a white dwarf 
as a result of thermonuclear burning of this matter. It is 
expected that during the explosion a mass of between a 
few 10~ 7 and a few 10 -5 M Q Q has been expelled with a 
velocity of 2760 — 3200 km s _1 . This mass interacts with 
the wind of RG providing conditions for acceleration of 
particles. It is estimated that such symbiotic Novae may 
appear within the Galaxy with the frequency of 0.5 — 5 
yr _1 0]. The binary system V407 Cygni is located at a 
distance of 2.7 kpc [l[ 

Recently, a 7-ray outburst has been observed from the 
direction of the binary system V407 Cygni by the Fermi- 
LAT telescope 0. The 7-ray emission was detected dur- 
ing the optical outburst, with the differential spectrum 
well described in the energy range 0.1 — 6 GeV by a power 
law with the spectral index —1.5 ±0.2 followed by the ex- 
ponential cutoff at E c = 2.2 ± 0.8 GeV. The emission was 
detected for ~ 15 days (during 2010 March 10-26) with 
the maximum at 3-4 days after the maximum of the op- 
tical emission [5[. Moreover, the 7-ray flux at the early 
stage (10 March - 14 March 2010) and late stage (14 
March - 29 March) was different by a factor ^2 while the 



spectral features remained unchanged within statistical 
limits [3]. V407 Cyg has been also observed by VERI- 
TAS at the later phase of the outburst during 2010 March 
19-26 @. Only the upper limit on the flux at the level of 
2.3 x 10~ 12 erg cm~ 2 s -1 has been reported at the energy 
of 1.6 TeV. 

The detected GeV 7-ray emission could have been pro- 
duced in the Inverse Compton (IC) scattering of the RG 
soft photons by rclativistic electrons or as a result of de- 
cay of pions produced in hadronic interactions 0. In 
fact, modeling of the Novae explosions shows that condi- 
tions are convenient for both radiation processes to occur 
. Both radiation processes can provide satisfactory fit- 
ting of the 7-ray spectrum measured by the Fermi Q. 
If GeV 7-ray emission is due to the hadronic processes, 
then neutrinos with ~ 10 GeV energies should be also 
produced as a result of decay of charged pions and muons. 
It has been argued that such low energy neutrinos (with 
energies ~ 0.1 — 100 GeV) from symbiotic Novae can be 
detected by the current and future experiments Q- 

In this paper we investigate whether symbiotic Novae 
can also produce TeV 7-rays and neutrinos detectable by 
the present and future instruments. We assume that the 
GeV emission observed by Fermi is due to the IC scatter- 
ing of the RG radiation by rclativistic electrons. Their 
maximum energies are determined by radiation energy 
losses. However, hadrons might be accelerated to signif- 
icantly larger energies due to the lack of saturation of 
the acceleration process by energy losses. These hadrons 
might produce TeV 7-rays and neutrinos in the interac- 
tion with the matter of the stellar wind. 



II. GENERAL SCENARIO 
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We apply a standard geometry of symbiotic binary sys- 
tem (sec Fig. [T|). We assume that both electrons and 
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FIG. 1: Schematic geometry of the symbiotic binary system 
V407 Cyg. The Nova shock wave propagates from the WD 
star (shock 1) and, as it propagates, starts to partially over- 
take the RG (shock 2). Particles (electrons, hadrons) are 
accelerated at the shock. An electron produces a 7-ray in IC 
scattering of a soft photon e (angle of interaction 0). Hadrons 
produce 7-rays and neutrinos in collisions with the matter 
from RG wind. 

hadrons can be accelerated at the Nova shock wave with 
comparable luminosities. In fact, this does not need to 
be true since in the case of the shock acceleration sce- 
nario the particle Larmor radius has to be larger than 
the thickness of the shock. This condition in general may 
not be fulfilled by electrons in the blast wave. Therefore, 
in principle electrons and hadrons might be accelerated 
at different sites and through different acceleration pro- 
cesses. For example, the magnetic reconncction, Fermi I 
and Fermi II processes can be considered here to operate 
at the shock wave or turbulent medium downstream of 
the shock. 

Taking into account the radiation energy density of RG 
and the density of matter of its wind (see Introduction) 
it is clear that the radiation process is expected to be 
most efficient in the vicinity of the RG. We assume that 
the observed GeV 7-ray emission is produced by elec- 
trons in the IC scattering of stellar radiation from the 
RG. In fact, a simple estimation of the optical depth for 
electrons on the IC process (in the Thomson regime) in 
the radiation field of RG close to its surface gives a value 
of the order of, tiq/t = "■rg-Rrg°t ~ 3.5, where tt-rq is 
the density of stellar photons determined by the surface 
temperature of this star T RG = 2.5 x 10 3 K= 0.25T 4 K, 
Rrg = 3.5 x 10 13 cm is the radius of RG, and ctt is the 
Thomson cross section. This indicates that electrons can 
lose significant amount of their energy through the pro- 
duction of 7-rays even on a distance scale comparable to 
the dimension of the star. In fact, electrons are captured 
in the vicinity of the star by the wind magnetic field, 
suffering efficient cooling during the whole lifetime of the 
flare ~15 days (i.e. corresponding to the propagation 
time of the shock through the separation of the stars in 



FIG. 2: Optical depths for 7-rays escaping from vicinity of 
the RG star. 7-rays are produced at the distance of 0.1 -Rrg 
(the dashed curves) or 0.5 -Rrg (the solid curves) from the 
surface and propagate along the radial direction away from 
the RG (the thin curves) or perpendicularly to it (the thick 
lines) . 

binary system. On the other hand, hadrons can be ac- 
celerated in the same region but to much larger energies 
since their energy loss process is not so efficient as in the 
case of electrons. These hadrons interact with the matter 
of the wind, losing a part of their energy through pion 
production (sec for estimates Q). 7-rays and neutrinos 
are produced as a result of decay of neutral and charged 
pions with energies which can extend up to TeV energy 
range. 

Since the TeV 7-rays are produced within the radiation 
field of the RG star, they in principle could be absorbed in 
the pair production process. However the optical depth 
for this process is much lower than 1 (see Fig. [2]), thus 
the loss of 7-rays due to absorption is negligible. 

Below we constrain the parameters determining the 
acceleration mechanism by comparing the 7-ray spectra 
produced by electrons in the IC process with the obser- 
vations of the GeV 7-ray emission from Nova V407 Cygni 
by Fermi-LAT. These constraints are farther used to de- 
termine the acceleration process of hadrons and predict 
the expected 7-ray and neutrino emission from their in- 
teraction with the matter of the stellar wind. 



A. Constraints on the acceleration of particles 

The gamma-ray spectrum observed from V407 Cygni 
by Fermi-LAT can be modeled by a power law with an 
exponential cut-off at energy E 1 ~ 2.2 GeV [||. If it is 
produced in the IC scattering of RG radiation, then in 
order to produce such 7-rays, electrons must have ener- 
gies 
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where m e is the electron rest energy, e = SksTaG = 2.6T4 
eV and /3 is the angle between the direction of electron 
and the thermal photon from the stellar surface. Close 
to the surface this angle is typically smaller than 0.57T. 

Let us consider the case of saturation of electron 
acceleration by the IC energy losses in the Thomson 
(T) regime occurring in the thermal radiation of RG. 
The cooling time scale of electrons on the IC pro- 
cess in the Thomson regime is estimated on i\c/t = 
3mg/4cr T cJ7 ra d-Ec ~ 5.4 x 10 4 i? 2 h /i? c s, where {7 rad is the 
energy density of RG radiation and the energy of elec- 
trons is expressed in GeV, and i? s h is in units of i?Rfi- 

The acceleration rate of electrons is parametrized by 
Pace = £cE/R L ~ l.l£-4-B GeV s"\ where R L is the 
Larmor radius of particles, the energy of electron E c is 
in GeV, the magnetic field B is in Gauss, £ = I0~ 4 £_4 
is the acceleration parameter, and c is the speed of light. 
By comparing the acceleration time scale, 

r a c cc = E c /P acc « lE e /(Z_ 4 B) s. (2) 

with the IC cooling time scale in the T regime, we esti- 
mate the maximum energies of accelerated electrons to 
be (see e.g. 0) 

Er x « 13(£_ 4 B) 1/2 ^ « 210(£_ 4 £0 1/2 i?sh GeV. (3) 

The value of the acceleration parameter determines the 
maximum energies of accelerated particles. In the case 
of the acceleration process occurring in the turbulent 
medium (Fermi II acceleration process) it can be esti- 
mated on £ ~ (uturb/c) 2 ~ 10 -5 , where Wturb ~ wn/3 = 
10 8 cm s _1 are the velocities of the turbulent irregulari- 
ties produced downstream of the shock, = 3 x 10 8 v$ 
cm s -1 is the velocity of the shock. We suppose that 
Fermi II acceleration process is more relevant for elec- 
trons which may not be able to fulfill the initial condi- 
tion for their injection in the shock acceleration scenario 
(Fermi I acceleration process). However, in the case of 
hadrons we assume that both mechanism may be rele- 
vant. 

Note that iJ™ ax should not depend on the distance 
from the star for the case of a radial structure of the 
magnetic field around RG. Such structure is expected in 
the case of a substantial stellar wind in the range of dis- 
tances characteristic for the separation of the stars in 
the V407 Cyg binary system. Therefore electrons can be 
injected soon after the Nova explosion. However the den- 
sity of radiation from the companion star, which creates 
a target for IC scattering, increases gradually when the 
shock approaches the star. We expect that the maximum 
production of 7-rays by electrons should be determined 
by the distance to the companion star and the velocity of 
the shock which might slowly decelerate lowering the ac- 
celeration efficiency of electron acceleration. As a result, 
the maximum 7-ray flux produced by electrons in the IC 
process is expected at some time after optical flash but 
before the crossing time of the binary separation by the 



shock, i.e. ~10 days. This is in agreement with the mea- 
surements of the evolution of the 7-ray flux measured by 
Fermi with the maximum ^3-4 days after optical flash. 

We expect that in the surroundings of the RG star 
the saturation of the acceleration process of electrons by 
the synchrotron energy losses is rather unlikely since it 
requires relatively strong magnetic fields in the stellar 
wind, > 40T 2 G, i.e above -B C rit ~ 2.5 G in the case of 
RG star. Such large values are not expected in the mag- 
netospheres of horizontal branch of the main sequence 
solar mass stars. The simple re-scaling of the surface 
magnetic field of the Sun (~3 G) to the dimension of the 
red giant phase (by a factor 500 -2 , assuming conserva- 
tion of the magnetic field flux, gives the value -~10~ 5 G 
and ~4x 10 -2 G for the active regions of its surface (the 
Suns spots with i? spo t ~ 10 4 G). The magnetic dynamo 
for red giants might be more efficient than the one oper- 
ating in the main sequence stars but it is rather unlikely 
that it will be able to produce a magnetic field strength 
that is larger by a few orders of magnitude. 

The comparison of the electron energy constrained by 
Eq. []] with that constrained by Eq. [3] allows us to esti- 
mate the product of the acceleration efficiency and the 
magnetic field strength at the shock, 

C-iB « 0.015(1 + cos /?)/i? 2 h G, (4) 

where i? s h is in units of Rrq. Note that this condi- 
tion is satisfied for the magnetic field strengths of the 
order of ~0.1 G for the required acceleration efficiency 
of electrons £ = 10 -5 . Such a magnetic field strength 
seems to be reasonable provided that the magnetic dy- 
namo operates more efficiently in the red giants than in 
the main sequence stars. Such value of the magnetic field 
is clearly below the critical value, -Bcrit, above which the 
synchrotron energy losses of electrons becomes more ef- 
ficient than their IC losses (note that most of the IC 
scattering will occur in the Thomson regime). There- 
fore, synchrotron energy losses can be safely neglected in 
our scenario. 

Below we assume that hadrons are also accelerated in 
the region of electron acceleration characterized by pa- 
rameters fulfilling the condition given by Eq. [?] 

B. Maximum energies of accelerated hadrons 

The acceleration time scale of protons, t^ cc , can be de- 
rived from an equation analogous to Eq. [2j Using the 
same value of £_4-B (given by Eq. [3} as for the acceler- 
ation of electrons (i.e. assuming that both populations 
of particles arc accelerated by the same processes in the 
same region), it can be estimated by, 

rP c = Ii? p /(£_4S)«70i? 2 h i?p/(l + cos/3) s, (5) 

where E p is the energy of proton expressed in GeV. The 
typical time scale of the flare (in case of no deceleration 
of the shock) is 

7t = (o - R sh ) x R R g/v n » 1.2 x 10 5 (a - R sh )/v 8 s,(6) 
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where a x Rug = 7i?RG = 15-5 AU is the separation 
between the stars. Tf is of the order of a week (similar to 
Fermi observations.) By comparing the acceleration time 
scale of the protons, given by Eq. [5j with the time scale 
of the flare, given by Eq. [6j we estimate the maximum 
energies of accelerated protons by, 

E^^1.7(a~R sh )(l + cos(3)/(Rl h v s ) TeV. (7) 

For the separation of stars in V407 Cygni given above and 
the shock arriving to RG surface (i.e. i? s h = 1), the max- 
imum energies of protons are estimated to be iJ™ ax ps 10 
TeV. These relativistic protons interact with the matter 
of the stellar wind producing charged and neutral pions. 
These products decay to 7-rays and neutrinos with en- 
ergies significantly larger than the 7-rays produced by 
the GeV electrons which scatter the stellar radiation (es- 
timated above). As shown above, the resulting 7-ray 
spectrum will not be strongly affected by the absorption 
in stellar radiation on the pair production process, even 
if it extends up tens of TeVs. 

We conclude that a two component 7-ray spectrum 
from the Nova V407 Cyg should be expected in such a 
model. The lower energy component (in the GeV en- 
ergy range) is produced in the IC of primary electrons. 
The second component (mainly TeV 7-rays) is produced 
by hadrons interacting with the wind of RG close to its 
surface. 



III. GAMMA-RAYS FROM THE NOVA V407 
CYGNI 

Applying the scenario considered above for the radia- 
tion processes, we calculate the IC spectrum produced by 
electrons which scatter RG radiation when the shock is 
close to the stellar surface. We allow electrons to cool 
completely in this radiation field. The energy distri- 
bution of electrons is selected such that the produced 
7-ray spectrum is consistent with the one observed by 
the Fermi-LAT. We consider two models for the electron 
spectrum: the monoenergetic spectrum (e.g. from mag- 
netic reconncction) and the power law with an exponen- 
tial cut-off (e.g. Fermi II process). The second case is 
defined as below. Since the cooling time scale of the elec- 
trons is clearly shorter than the time scale of the flare (as- 
sumed equal to the propagation time of the shock to the 
RG), we estimate that the injected differential spectrum 
of electrons should have the form, dN e /dE c oc E~ x - 5 , in 
order to produce the 7-rays in the IC process with the 
spectral shape cx -E" 1 ' 75 . From the normalization of the 
calculated 7-ray spectrum to the observed spectrum, we 
estimate the total power in the electron spectrum. 

We also assume that similar spectral index and to- 
tal power is characteristic for the spectrum of relativis- 
tic protons since they are expected to be accelerated si- 
multaneously by this same acceleration mechanism (note 
however that they are accelerated to different maximum 
energies, as their energy loss processes are different). The 



collision time scale of relativistic protons with the matter 
of the stellar wind (for p — p — > tt process) is, 

Tp_ p = (n w c(Tp_p) _1 w lO 7 ?^ 1 s, (8) 

where the density of matter in the wind is n = 10 8 ns 
cm and the cross section for pion production in colli- 
sions of protons in matter is cr p _ p w 3 x 10~ 26 cm 2 . The 
density of RG wind, close to the stellar surface, is esti- 
mated from n = M / {An R\ G m v v KG ) «7.8 x 10 8 A/_ 6 . 5 /V6 
cm -3 , where m p is the mass of proton, the mass loss rate 
of the Red Giant is M =3 • KT 7 A/_ 6 . 5 M yr _1 and the 
velocity of the wind from RG is tjrq =10 6 Vg cm s _1 . 
We conclude that a significant amount of energy (up to 
a few tens %) can be lost by relativistic hadrons during 
propagation of the shock through the binary system. 

We calculate the 7-ray spectra from IC scattering of 
the thermal photons from RG with the help of a dedi- 
cated Monte Carlo simulations code. Since the electrons 
with energies of a few tens of GeV (close to the cut-off of 
the spectrum) will no longer interact in the pure Thom- 
son regime, we use the full Klein-Nishina cross-section 
for the IC scattering in the simulations. On the other 
hand the interactions of the relativistic protons with the 
stellar wind matter are studied with the Geant 4 sim- 
ulation toolkit [{J [l(|. We apply the QGS model of 
hadronic interactions which allows us to perform calcula- 
tions for protons with energies up to 100 TeV. While mul- 
tiple hadronic interactions of a single proton are rather 
rare in the above described scenario (of the order a few 
%), they are still automatically taken into account during 
the simulations. 

As a first example, we consider an acceleration of mo- 
noenergetic electrons with energy ~ 25 GeV at the dis- 
tance 0.5 -Rrg from the surface of the RG star. Simul- 
taneously protons with energy ~ 10 TeV are accelerated 
at the same shock as predicted by Eq. 7. The electrons 
will cool down completely in the radiation field of the 
RG producing GeV 7-rays. We normalize the electron 
spectrum to reproduce the total 7-ray flux measured by 
the Fermi-LAT. We assume that the total luminosity of 
the accelerated protons is equal to that one of electrons. 
We investigate the 7-ray spectrum obtained from the de- 
cay of 7T° from hadronic interactions. Moreover, the e 
created in the decay of can also scatter the stellar 
radiation producing additional component in the 7-ray 
spectrum. Those three components are shown and com- 
pared with sensitivities of present and future Cherenkov 
telescopes in top panel of Fig. [3] 

The luminosity of electrons, L e , required to explain the 
GeV emission observed by Fermi-LAT is ~ 4.2 x 10 41 erg 
(assuming they are accelerated over the 14 day period), 
consistent with value ~ 4 x 10 41 erg obtained in modelling 
of IC emission in [i| . With a plausible observation time 
of 20h, TeV 7-ray emission from the hadronic component 
could be detected from a comparable Nova outburst by 
present Cherenkov telescopes so long as the luminosity of 
accelerated protons, L p , is above ~ 0.3 of the accelerated 
electron luminosity. With the future CTA instrument 
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FIG. 3: The comparison of the Fermi-LAT 7-ray aver- 
age spectral energy distribution (SED) over the period 2010 
March 10-29 (black points and upper limits, Q) with the cal- 
culations of the IC 7-ray spectrum produced by electrons and 
the TeV 7-ray spectrum produced in hadronic collisions. Mo- 
noenergetic electrons with energy 25 GeV and protons with 
energy 10 TeV (top panel) or electrons following a spectrum 
dN/dEe oc E~ 1,5 exp(_B e /20 GeV) and protons with spectrum 
dN/dE p oc Ep 15 exp(E p /8TeV) (bottom panel) are acceler- 
ated in turbulent medium. Those particles are injected at 
the distance 0.5 Rrq from the surface of RG. The thick black 
line is the total predicted 7-ray emission composed of three 
components: IC scattering of the injected electrons (the thin 
solid line), n° decay (the thin dashed line) and IC scattering 
of the electrons created in decay (the dotted line) Thick 
gray lines show the differential sensitivity (scaled down to 20 
hrs) of the MAGIC telescopes (the solid line, pj]]) and fu- 
ture CTA observatory (dashed line, Q2])- The empty square 
marker with an arrow is the flux upper limit from VERITAS 
observations during the period 2010 March 19-260]. 



this component can be detected even if L p > 0.03L e . The 
upper limit on the 7-ray flux from V407 Cygni derived by 
the VERITAS telescopes [|| is below the predicted 7-ray 
flux for L p = L e , which is based on the normalization 
to the average GeV flux observed by Fermi-LAT. Simple 
comparison suggests L p < 0AL e . Note however, that the 
observations by VERITAS were performed 5 days after 



the peak of the GeV emission. Thus, direct comparison 
depends on the exact time evolution of 7-ray emission 
produced by electrons and hadrons. A simple investiga- 
tion of the GeV 7-ray light curve (Fig. 1 in Q) indicates 
that the power in relativistic hadrons should be not larger 
than the power in relativistic electrons assuming similar 
time evolution of the 7-ray emission in these two energy 
ranges. 

Secondary pairs, produced in hadronic interactions, 
have some small effect on the overall spectrum. The IC 
scattering of these pairs softens the total 7-ray spec- 
trum in the sub- TeV range. It also provides a bridge 
emission between the low and high energy components. 

As mentioned above, we also consider a scenario in 
which both, the electrons and the protons, follow a 
power-law energy distribution with the same spectral in- 
dex but with different exponential energy cut-offs. We 
assume a spectral index of 1.5 and cut-offs determined 
by Eq. [2] and Eq. [71 The total energy in accelerated 
electrons ~ 1.5 x 10 42 erg (accelerated over the 14 day 
period) needed to explain the GeV emission is slightly 
higher than in the injection of monocncrgctic electrons. 
As before we assume that the total luminosities of pro- 
tons and electrons are comparable. The Fermi-LAT spec- 
trum can be fitted for energies above ~ 0.3 GeV, however 
the predicted flux below this energy is overestimated (see 
bottom panel of Fig. [3]). This may be e.g. due to higher 
minimal energy of injected electrons, or only partial cool- 
ing of the electrons. Please note however that measure- 
ments in this energy region have high both statistical and 
systematic errors. As the 7-ray spectrum in sub- TeV en- 
ergy range is softer in this case than for monoenergetic 
protons, the effect of the IC scattering of e ± produced 
in the hadronic interactions is less pronounced. The ob- 
servations of VERITAS yield a similar limit on the to- 
tal luminosity of accelerated protons as in the case of 
monoenergetic injection. However, since the VERITAS 
observations were performed with a larger than normal 
energy threshold (due to high zenith angle of the source 
during the observations), they did not allow to explore 
the possible sub- TeV emission, which is predicted to be 
strong in the power-law injection scenario. 



IV. NEUTRINOS FROM V407 CYGNI 

We consider two models for acceleration of hadrons 
differing in acceleration efficiency. In the first one (de- 
scribed in detail in the previous section), both hadrons 
and electrons are accelerated either in magnetic reconnec- 
tion process or the Fermi II acceleration process occur- 
ring in the turbulent medium downstream of the shock. 
In the second model, we assume that hadrons are ac- 
celerated on the shock wave via the Fermi I process. 
In this case the acceleration efficiency for hadrons (de- 
scribed by the value of £) is expected to be significantly 
larger than that for electrons estimated above for the tur- 
bulent medium. In the Fermi I process £ can be as large 
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FIG. 4: v + v emission from hadronic interaction of 
lOTeV protons (thin lines) or with spectrum dN/dE p oc 
Ep 1 ' 5 exp(B p /8TeV) (thick lines). Black dashed line are the 
muon neutrinos, and black dotted line are the electron neutri- 
nos. For comparison we show the atmospheric muon neutrino 
flux in a viewcone of 1° radius with a thick solid gray line 
[l3| . and the sensitivity of the IceCube detector (scaled to 
14 day observation time) with a dashed gray line (IceCube 
design report). 



as ~ v-m/c. 

For both proton injection models considered above, we 
also calculate the neutrino spectra produced in the de- 
cay of 7r ± (and subsequent fi ± ). In Fig. |H we show the 
predicted electron and muon neutrino v + v fluxes for 
the monoenergetic and the power law injection of TeV 
protons (see the first example in Section llllj) . They are 
compared with the atmospheric neutrino background and 
the sensitivity of the IceCube telescope. Note that for the 
considered model the spectra are clearly below these lim- 
its for the case L p — L c . As we mentioned above, L p is 
not likely above L due to the TeV 7-ray observations 
reported by VERITAS @. Therefore, we conclude that 
detection of a single outburst by neutrino instruments of 
the IceCube type is not very likely. 

As the second model we consider the case of acceler- 
ation of protons with the power law spectrum up to an 
order of magnitude higher energies than in the previous 
case and with the luminosity being 10 times larger than 
that of accelerated electrons. Such a scenario might be 
consistent with the upper limits on the TeV 7-ray flux 
derived by VERITAS from observations of V407 Cygni, 
provided that the time dependence of the acceleration 
of electrons in the turbulent medium (as considered in 
the first scenario) and protons on the shock wave in the 
Fermi I type mechanism also differs significantly. For ex- 
ample, this might be the case when the magnetic field at 
the blast wave is dominated by the magnetic field from 
the White Dwarf. Then, the maximum energies of pro- 
tons are expected to be larger when the magnetic field 
at the shock is stronger, i.e. at the very early phase of 
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FIG. 5: SEDs of 7-rays (the upper panel, like in bottom panel 
of Fig. [3| and v + v (the lower panel, like in Fig. 4) but for 
the case of efficient hadronic acceleration in the first order 
Fermi acceleration process. Injection of electrons dN/dE oc 
Eg exp(_E e /20 GeV) and protons with spectrum dN/dE oc 
Ep 1 ' 5 exp(_E p /80 TeV) occurs at the distance 0.5 -Rrg from 
the surface of RG. The total power in relativistic protons is 
10 times larger than that one in electrons. 



Nova explosion. In the latter phase, when the magnetic 
field at the shock drops, protons are not accelerated to 
such large energies due to their escape from the shock re- 
gion. On the other hand, electrons start to be efficiently 
accelerated when the downstream region of the shock be- 
comes more turbulent, i.e. at the phase of the explosion 
corresponding to the maximum 7-ray emission observed 
by Fermi-LAT. 

In Fig. [5] we show the expected 7-ray and neutrino 
spectra for the case of the power-law acceleration of elec- 
trons and protons. In this scenario the neutrino spectra 
overtake the atmospheric background in the multi-TeV 
range by an order of magnitude. Moreover such a neu- 
trino outburst is above the sensitivity limit of the Ice- 
Cube telescope. However, in this case strong 7-ray emis- 
sion from pion decay should be easily detected by the 
current IACTs even in relatively short time scale. The 
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VERITAS upper limit at 1.6 TeV excludes such a possi- 
bility during the latter phase of the 7-ray outburst (i.e. 
after March 19th 2010). Therefore, in order to be consis- 
tent with those observations hadrons should be efficiently 
accelerated only during the early phase of the Nova out- 
burst. 



V. CONCLUSION 

We have considered a two-component scenario for the 
high energy processes during symbiotic Nova explosions, 
proposing that GeV 7-ray emission is produced by elec- 
trons and the accompanying TeV 7-ray and neutrino 
emission is produced by hadrons. The GcV 7-ray spec- 
tral shape seems to favor monoencrgetic electron injec- 
tion, such as can be obtained e.g. in magnetic reconnec- 
tion process, rather than power-law electron production 
typical in e.g. Fermi II acceleration. Two models are con- 
sidered for the acceleration of the hadrons. In the first 
model, hadron acceleration process is tuned to the accel- 
eration process of electrons which is constrained by the 
GeV 7-ray observations. This model predicts relatively 
low fluxes of TeV 7-rays which are marginally consistent 



with the recent upper limits by VERITAS observations 
provided that similar power goes on the acceleration of 
electrons and protons. With future CTA project the TeV 
component in gamma-ray spectrum can be detected even 
if the luminosity of protons is as low as 3% of the one of 
electrons. The expected neutrino fluxes are below sensi- 
tivities of the IceCube telescope. In the second model, 
hadrons are accelerated independently from electrons in 
a more efficient mechanism (Fermi I for protons in con- 
trast to Fermi II or reconnection for electrons). This 
rather optimistic scenario predicts strong fluxes of TeV 
7-rays which in fact may only appear during the early 
stage of the Nova flare (before March 19th) due to the 
constraints set by the VERITAS upper limits. Also the 
accompanying neutrino burst at TeV energies might be 
detectable in this case by the IceCube telescope. 
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